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Abstract 

We consider a particle which is randomly accelerated by Gaussian white noise 
on the line < x < 1, with absorbing boundaries at x = 0, 1. Denoting the 
initial position and velocity of the particle by xq and vq and solving a Fokker- 
Planck type equation, we derive the exact probabilities qo(xo, vq), qi(xo,Vo) 
of absorption at x = 0, 1, respectively. The results are in excellent agreement 
with computer simulations. 

Permanent address: Department of Physics, Temple University, Philadelphia, PA 19122. 
PACS 05.10.Gg, 02.50.Ey, 05.40.-a 



1 



A well known topic in random walk theory |l[ is the problem of the "gambler's ruin." 
Initially the gambler has an amount of money Xq and the bank the amount 1 — Xq. The 
gambler flips a coin repeatedly, randomly winning or losing the increment e. The game ends 
when the gambler's funds reach or 1. The problem is to compute the probability qo{xo) 
that the gambler loses everything. 

The problem is easily solved. Since qo(xo) = k[Qo( x o + e ) + %( x o ~ e )]? 

d 2 q (x ) 



(1) 



in the limit e — > 0. /,From equation ([!]) and the boundary conditions qo(0) = 1, ?o(l) = 0, 

q (x ) = l—x . (2) 

As the starting capital increases from to 1, the probability of the gambler's ruin decreases 
from 1 to 0. 

Instead of the gambling scenario one could equally well imagine a particle making a 
random walk with infinitesimal steps ±e on the x axis, with initial position < xq < 1. In 
the course of time the particle eventually arrives at x — or x — 1. The quantities qo(xo) 
in equation (Q) and gi(x ) = 1 — ^(^o) represent the probabilities that the particle first 
reaches the edge of the interval at x = and x — 1, respectively. Alternatively, we could 
impose absorbing boundary conditions and interprete go^o) an d qi(%o) as the probabilities 
of absorption at x = and at x — 1. 

In this Letter we also consider a particle on the finite interval < x < 1, but we assume 
that the changes in the velocity rather than the position of the particle are random. The 
particle moves according to the Langevin equation 

^ = V(t), (3) 
where the acceleration r](t) has the form of Gaussian white noise, with 

(17(f)) = 0, (v(t 1 )v(t 2 )) = 2S{t 1 -t 2 ). (4) 
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Imposing absorbing boundary conditions, we derive the probabilities qo(xo, vq), qi(xo,vo) 
of absorption at x = and at x — 1, respectively, as functions of the initial position and 
velocity 

The quantity qo(xo,vo) can also be interpreted as the probability of a gambler's ruin, 
but the game is different. The gambler has an amount of money x(t) at time t and the 
bank the amount 1 — x(t). Money is transferred from the bank to the gambler at a rate 
v = dx/dt, which may be positive or negative. At regular infinitesimal intervals the gambler 
flips a coin, randomly increasing or decreasing the rate v by the increment A. The game 
ends when x reaches or 1. The quantity qo{xo,vo) is the probability that a gambler with 
initial conditions x ,v loses everything. 

In the case of a random walk on the x axis, the probability density P(x, Xq, t) at time t of a 
particle which is initially at Xq obeys the diffusion equation. For a particle which is randomly 
accelerated according to equations (|3|) and (P, the probability density P(x, v; xq, v q; t) in 
the phase space (x, v) satisfies the Fokker-Planck equation |2[ 

( d d d 2 \ 

\¥t +l Tx-^) P ^ V ' M]t)= ^ (5) 
corresponding to diffusion of the velocity, with initial condition 

P(x,v;x ,v ; 0) = 5(x - x )S(v - v ). (6) 

In analogy with the discussion leading to differential equation ([I]) for q (xo), let us 
consider a discrete dynamics in which the velocity v changes by ±A with equal probability 
at time intervals r. For this dynamics 

P(x, v; x , v ; t) = ~[P{x -vr,v + A; x , v ] t - r) + P(x -vt,v- A; x , v ; t-r)], (7) 
qo(x ,v ) = ^[q (x + v t,v + A) + q Q (x + v t,v - A)]. (8) 

Expanding equations (0) and (^j) in r and A, dividing by r, and taking the limit r = |A 2 — > 
gives us a "poor man's" derivation of the Fokker-Planck equation (Q) and the corresponding 
differential equation 
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for the probability of absorption at x = 0. 

To solve equation (pi) with the absorbing boundary condition 



g (0,vo) = l, ^o<0, (10) 



and the requirements 



qo(x ,v ) = qi(l- x ,-v ), (11) 
<lo(xo, v ) + qi(x , v ) = 1, (12) 

of reflection symmetry and total probability equal to 1, we first make the substitution 

ip{x,v) = q {x, -v) - -. (13) 



Expressed in terms of ip(x,v), equations ©-(12) take the form 

( d d 2 ^ 



dx dv 2 



i/;(x,v) = 0, (14) 



= ^>0, (15) 
i/j(x, v) = — if)(l — x, — v). (16) 

Masoliver and Porra [§] have shown how certain Fokker-Planck type equations on the 
finite interval < x < 1 can be solved exactly. They derived an exact result for the 
average time T(xq,Vq) a randomly accelerated particle with initial conditions Xq,Vq takes 
to reach a boundary of the interval. The probability that the particle has not yet reached 
a boundary after a time t decays as e~ Et , as discussed by Burkhardt He obtained E 
numerically with an approach similar to and related it to the confinement free energy of a 
semiflexible polymer in a tube. In another application inspired by [[|, Burkhardt, Franklin, 
and Gawronski || calculated the equilibrium distribution function P(x, v) of a randomly 
accelerated particle on the line < x < 1 undergoing inelastic collisions at the boundaries 



The function ip(x,v) satisfies the same steady-state Fokker-Planck equation fll4]) as the 
quantity P(x, v) considered in || and has the same Green's function solution 



ip(x, v) 



3x Jo 
1 



duu^ 2 e-^ +u ^ 9x J_ 1/3 f L " n " | V(0,«) 



V 9x 

-«3/90r-y) 



v > 



(17) 



3V3r(f) Jo y (x-y) 2 / 3 dv 

derived in ||. Equation QT7D only holds for positive v. For negative v, if)(x,v) can be 
obtained from equation flTT| ) using the antisymmetry ( |TE| ) under reflection. 

Equation ( 17) determines ip(x,v) for all x > and f > from ip(0,v) and dip(x,0)/dv. 
The first of these functions is given in equation (|15f) . To determine the second, we set t> = 



in equation (|17D , which yields 
1 



^(x, 0) 



3V3r(? 



x" 2 / 3 / 4 U e-"^(0, U )- , 

jo Jo (x — yY' 6 ov 



dy dip(y, 0) 



(18) 



Then, substituting equation (|18|) in the relation i[>(x, 0) +'0(1 —x,0) = 0, which follows from 
(PH), and using dip(y, 0)/dv = dip(l — y, 0)/dv, also a consequence of ([16]) , we obtain 

dy dtp{y,0) 



\x — y\ 2 ^ dv 



duu 



e ~u a /9x e -u 3 /9(l-x) 



+ 



x 2/3 (l- X )2/3 



^(o,«). 



(19) 



The solution to integral equation (]TD[) , derived, following [0], in Appendix B of |J, is given 
by 

dip(x, 0) 



<9t> 



dit ti [R(x, u) + i?(l — x, u)\ ip(0, u) 



(20) 



where 



R(x, u) 



u l/2 e -u i /9x 



H F (_l 5 

I _ x U/6 11 V 6' 6' 9: 



-x) 



9x J ' 



(21) 



3 5 / 6 r(i)r(|) x 7 i Q (i- 

and i -Fx (a; 6; z) is the confluent hypergeometric function @,|[. 

Equations (|TTD , (p0|) , and (f2~l|) determine ip(x,v) for all x and v from -0(0, v) for t> > 0, 
which is known from the absorbing boundary condition (|1|). Substituting equations ( |15|) 
and fl2T|) in ([21]) leads to 
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and from (|Tj|,(|17|), and (||) 

= i- *L *MW |s(1 -» r " (23) 

Rewriting equation (|23|) in terms of qo(xo,Vo) using (|TT| ) -([T3| ) , we obtain our main result 

1 r 1 e - v o/ 9 (y~ x o) 

q (x ,v ) = l-q {l-x ,-v ) = — dy — [y(l - y)]~ 1/6 , v > 0, (24) 

in jxq [y — xq) 1 

analogous to the solution (|2|) of the traditional gambler's ruin problem. 

For Xq — 1 equation ( p4|) reproduces the expected result <7o(l,i>o) = 1 — 9o(0, — vq) = 0, 

fo > 0, corresponding to the immediate absorption of a particle that is initially at either 

boundary with velocity directed outward from the interval < x < 1. For xq = and Vq = 

the integral in equation can be evaluated, yielding 

2'3 2 / 3 1/2 771/1 7 1 3\ 

YfTV^O lulls' 6? -9 V 0)i 



q (0, v ) = l- q (l, -v ) = 1 - — j- v 1 ^ f ; -§i; 3 ), t; > 0, (25) 

q (x , 0) = 1 - q (l - x, 0) = 1 - ^1 X V 6 ^(i, | ; I ; Xq ). ( 26 ) 



^6' 

Here i-F\(a; b; z) and 2 c; 2) are the confluent and ordinary hypergeometric functions 

II- 

The probability (70(^0, u o) of absorption at the origin, obtained from equation (^4j) by 
numerical integration, is shown in figure 1. The probability decreases monotonically as xq 
increases with fixed v and as v increases at fixed x , as expected. The quantity qo(x , v o) 
is a nonsingular function of (xq,Vq) except at the two boundary points (0,0) and (1,0). 
The curves for xq = 0.0,0.1,0.3,0.5 become smoother near v = as Xq increases, and for 
x = 0.5, qo(x ,v ) — I is an odd function of v , as implied by equations ([□]) and (P"2"|). 



The points in figure 1 show the results of computer simulations, which clearly are in 
excellent agreement with the analytical results. Our simulation routine is based on the 
exact solution |TIJ 



r, / f 3 

P bee {x,v;x ,v ]t) = —7, expj- — 



(X-X - V t)(x - Xq - Vt) + -(v - V )H 2 



(27) 



of the Fokker-Planck equation flj^) with initial condition (|6]) in the absence of boundaries. 
Trajectories with the probability distribution Pf ree (x n+ i, v n+ \, x n , v n ; A n+1 ) given by fl2"T|) are 
generated using the algorithm 

/A 3 \ 1/2 

x n +i = x n + v n A n+1 + ( — j (^ Sn+1 + v 7 ^ r n+1 ) , (28) 
^n+i = + (2A n+1 ) 1/2 r n+1 , (29) 

where x n and t> n are the position and velocity of the particle at time t n , and A n+1 = t n+ i—t n . 
The quantities r n and s n are independent Gaussian random numbers such that 

(r n ) = (s n )=0, (r 2 n ) = (4) = 1. (30) 

In the absence of boundaries there is no time-step error in the algorithm, i.e., the A„ may 
be chosen arbitrarily. Close to boundaries small time steps are needed. 

To derive a quantitative criterion for an acceptable time step, we begin with the averages 

(x(t)) =x + v t, ([x(t) - (x{t))] 2 ) = h\ (31) 



implied by the distribution function (pT|) . At time t the particle coordinate x has a Gaussian 
distribution, with a maximum at x = xq + vot and the root-mean-square width ft d . 



The effect of the boundaries on the propagation is negligible if the Gaussian peak lies almost 
entirely within the interval < x < 1. This is certainly the case if, say, 

< x + v t ± 5t 3/2 < 1 (32) 

Over the range of velocities encountered in our simulations, any t which satisfies the simpler, 
more stringent condition 

t<—x (l-x Q ) (33) 

also satisfies (|32"D . 

Keeping inequality (|33D in mind, we performed our simulations with the time step 

A n+1 = 1(T 5 + 10" 1 a; n (l -x n ). (34) 



The time step decreases as the particle approaches the boundary and has the minimum value 
10 -5 . It is necessary to have a small nonzero minimum value. Otherwise the particle never 
arrives at the boundaries. Our results for the absorption probability qo{xo,vo) are averages 
based on 10 5 trajectories for each set of initial conditions xq,Vq. 

Finally we note that qo(xo,Vo) may be derived from another general Green's function 
solution Q of the Fokker-Planck equation ([14]), 



v 1 / 2 r°° , , , ,, , /2?; 3 / 2 w 3 / 2 \ 
*{x,v) = V —J o duu^e-^^ J 1/3 [^—) mu) 

v ex e -v 3 /9(x-y) 



+ wmjo dy v^w i){v ^ v>0} (35) 



different from fllTD - By substituting equation (R3) in (II), one can check that the Fokker- 
Planck equation is indeed satisfied. On the lines x = and v — equation (R3) reduces to 
the identities ^(0,^) = ip(0,v) and ip(x, 0) = ?/>(:e,0), respectively. 

Equation (|35| ) determines i/j(x,v) for all x > and v > from i[)(0,v),v > and 



if)(x, 0). The first of these functions is given in equation (|T^). To determine the second, we 
differentiate equation (|35f) with respect to v and then set t> = 0, which yields 

dip(x,0) _ 1 

<% ~ 3 2 / 3 r(|) x 

fx rlii f)ih(n f\\ 

(36) 



3 

oo 







/o (x — y) 1 / 3 dy 

For the absorbing boundary condition ( |T5| ) the first two terms on the right-hand side of ( |36| ) 
cancel. Substituting equation fl3"6j ) in the relation dip(x, 0)/dv — dtp(l — x,0)/dv = 0, which 
follows from ([0]), using the invariance of di/j(y, 0)/dy under y — > 1 — y, and integrating with 
respect to x yields 



Cdy 
Jo 



I |2/3^(^ ) , / QV \ 

~~ — d = co ^ ' 



1 This solution may be derived by slightly modifying the derivation in Appendix A of |||. Setting 
v = in equation (A3) of ||, solving for W(s), and reinserting the result in (A3) with v ^ yields 
the Laplace transform of the new solution (|35|) . 

8 



The function ip(x,0) given in equations (p~3|) and fl26| ) satisfies equation (|37|). Substituting 
this ^(a;, 0) and -ip(0, v ) = \ into equation (p5|) , integrating, and using (|TTD-(|T3D , we obtain 

q (x , v ) = 1 - q (l - x , -v Q ) 

=T<rl d!/ (^)^ (1 - !/) 2Fi( ^- ;1 - y) ' "° >a (38) 

With the help of the identity 



for arbitrary f(y), one can convert expression (38) for qo(xo, Vq) to the simpler form (p4|). 

The second of the two Green's function solutions (|17|) , (|35l) looks simpler than the first, 
since no derivatives of ip appear on the right-hand side, but our main result fl24|) for go(^o, ^o) 
is obtained more easily from (flTl). 
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FIGURES 
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FIG. 1. The probability qo(%o, v o) of absorption at the origin. The solid lines show the exact 
result given in equation (p4|). The points are the results of our computer simulations. The data 
points have a statistical uncertainty ±5qo with \5qo\ ~ 0.001. 
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